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Starting from the recognized fact that oscillations of moments with rank and shoulder structure in the multiplicity 
distribution have the same origin in the full sample of events in e + e~ annihilation, we push our investigation to 
the 2-jet sample level, and argue in favor of the use of the negative binomial multiplicity distribution as the 
building block of multiparticle production in e + e~ annihilation events. It will be shown that this approach leads 
to definite predictions for the correlation structure, e.g., that correlations are flavour independent. 



1. FRAMEWORK AND DEFINITIONS 

This work is a step in the direction of the de- 
scription of multiplicity distributions (MD's) and 
correlations functions within a common frame- 
work: this is the approach we chose in order to 
investigate the dynamical mechanism of multipar- 
ticle production. In particular, in this work we 
explore the relationship between the multiplicity 
distribution, P n , i.e., the probability of producing 
n charged final particles, and the g-particle cor- 
relation function C q {y\, . . . , y q ), using e + e~ an- 
nihilation data. 

We will use the following moments of the MD: 
a) the factorial moments: 



F q = J2n(n -!)■■■ (n-q+l)P n 



b) the factorial cumulant moments: 



(1) 



9-1 



i=l 

c) their ratio H, 

K„ 
q F 



(2) 
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and will study these moments as a function of 
the order q (see for more details). While the 
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above definitions are valid in any restricted region 
of phase space, the work reported here refers only 
to full phase space. 

It should be noticed that the following relation 
holds: 



K q = J dyi 



■dy q C q (yi, ...,y q ) 



(4) 



where the integration is over the full phase space, 
and C q (y\, . . . , y q ) is the correlation function ob- 
tained from the inclusive m-particle densities 
(to = 1, . . . , q) by cluster expansion 

The discussion in the following sections is based 
on a detailed analysis of published experimental 
data. A few points should be made clear in ad- 
vance: 

a) Because of charge conservation, the number 
of particles in full phase space must be even, so 
all MD's in the following must be understood as 
'in their even component only', e.g., 



should be read as 



P 



(NBD) 



(/i, k) if n is even 
if n is odd 



(5) 



(6) 



where a proportionality constant serves the pur- 
pose of normalization. This procedure is normally 
used in the literature (|]. 

b) Because the experimental data sample is fi- 
nite, the tail of the distribution cannot be fully 
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Figure 1. Final charged particles MD for the full 
sample of events at y/s = 91 GeV from the Del- 



phi Collaboration (diamonds) and fits with eq. 
(dotted line, open circles in the residuals) and 
(solid line, filled circles in the residuals). 



Figure 2. H q ratio from the MD's in figure [j] 



than 'solid results'. As predictions, they are in- 
deed quite intriguing and our approach should 
therefore be considered as a strong suggestion to 
the experimental groups. 

d) For the same reason discussed in c), the er- 
ror bars we plot here for the ratio H q should be 
considered indicative: the method of varying ran- 
domly each P n within its standard error with a 
Gaussian distribution many times in order to cal- 
culate the variance of the H q should be applied to 
the uncorrelated P n (i.e., before unfolding); how- 
ever, even if we apply it to published data, it 
agrees with experimental results H . 



sampled, and the published MD's are indeed trun- 
cated at some value; all analytical MD's in the 
following are understood to be truncated at the 
same value as the corresponding experimental 
ones. The importance of truncation on the anal- 
ysis of the ratio H q is well known and has been 
investigated in ||. 

c) Fits performed on the published MD's suf- 
fer from the fact that the matrix used to go from 
the measured MD to the reconstructed MD in the 
unfolding procedure (see, e.g., Q) is not usually 
published. Thus a correct, unbiased regression 
cannot be made; the fits we describe here should 
then be seen more as 'solid predictions' rather 



2. THE FULL SAMPLE OF EVENTS 

We begin by examining the full sample of 
events at y/s — 91 GeV. Data on MD's have been 
published by all Collaborations, but in order to be 
consistent and concise we will only use the data 
published by the Delphi Collaboration j| (ref- 
erence |j| contains a more detailed analysis). In 
||, the negative binomial distribution (NBD) is 
fitted to the data; it is worth recalling the explicit 
expression for the NBD: 

Pi NBD) (M,fc) = 

k(k + l)'--(k + n-l) fi n k k 

n\ (/i + k) n + k ( ' 
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where the two parameters n and k are related to 
the first two moments as 



and k 



(8) 



The description of the data by the NBD is clearly 
unsatisfactory from the point of view of the chi- 
square (see table 0a), and also from the point of 
view of the residuals (shown in figure |l| with open 
dots): there is a clear structure which relates to 
what has been called 'shoulder structure' in P n 
Q. The corresponding situation with the ratio 
H q is shown in figure |^, from which it is evident 
that a single truncated NBD (dotted line) cannot 
describe the oscillations of H q with q. 

It was shown by the Delphi Collaboration j(| 
that the shoulder structure can be explained in 
terms of the superposition of the MD's resulting 
from classifying the events according to the num- 
ber of jets, and that the MD in each class is well 
fitted by a NBD. Moreover, the goodness of this 
description is not influenced by varying the pa- 
rameter of the jet-finding algorithm. 

It is then natural to try a description of the 
full sample of events with the weighted sum of 
two NBD's ||, one corresponding to 2-jet events 
and the other to 3-jet events (4-jet events' con- 
tribution is negligible). This is implemented by 
forcing the weight to be equal to the fraction of 2- 
jet events as experimentally measured in [|| . The 
resulting MD is thus 

Pn = ^-jet^n" 8 (/^2-jet) &2-jct) + 

(l-a 2 . jet )Pi NBD )(M3-jet,fc3-jet) (9) 

This fit is shown with a solid line in figure |l| for 
c*2-jet = 0.767, corresponding to y m i n = 0.06 in 
the jet-finder JADE; the parameters of the fit are 
shown in table |l|b. In the same figure we also no- 
tice that the structure of the residuals has greatly 
improved, and in figure || that the H q ratio is also 
well described (solid line). Furthermore, the pa- 
rameters we find from our fit are consistent with 
those obtained experimentally fitting the 2- and 
3-jet samples separately with single NBD's. Anal- 
ogous result can also be obtained analyzing the 
SLD and OPAL data §; L3 data @ are consis- 
tent with SLD's. 



Table 1 

Parameters of the MD's fitted to the data as 
shown in the figures. Part a) refers to figure |l| 
(dotted line); part b) to figure [l] (solid line); part 
c) to figure H (dotted line); part d) to figure || 
(solid line). 



a) 




21.4 ±0.9 




k 


24.3 ±0.7 




X 2 /NDF 


80/34 


b) 


M2-jct 


19.4 ±0.2 




&2-jet 


52 ±6 




M3-jet 


27.3 ±0.3 




&3-jet 


61 ± 10 




X 2 /NDF 


12/21 


c) 


M 


18.5 ±0.1 




k 


57 ±3 




X 2 /NDF 


20/16 


d) 


Mlight 


17.2 ±0.2 




Mhcavy 


22.0 ± 1.6 




k 


145 ± 53 




X 2 /NDF 


13/16 



3. THE 2-JET SAMPLE OF EVENTS 

It is a striking result, and a lesson to be learned 
from the analysis mentioned in the previous sec- 
tion, that one can resolve puzzling structures by 
investigating them on a deeper level. It is a pre- 
diction of the same analysis that the 2-jet sample 
should be well described (in all its aspects) by 
one NBD. However, this turns out not to be en- 
tirely true, as shown in figure |3|, which reproduces 
the 2-jet data and (single) NBD fit by Delphi [f|. 
While the chi-square is satisfactory (see table |l z) , 
the residuals (open dots in the figure) show some 
structure, although much less pronounced than 
in the full sample; the ratio H q (figure ||) still 
shows oscillations, albeit an order of magnitude 
less deep than in the previous case. 

A few results coming from experiments will 
help us investigate deeply the structure of 2-jet 
events. First, the result by OPAL || on forward- 
backward correlations: they are weak, and in the 
bulk they are explained by the superposition of 
events with a fixed number of jets; the residual 
correlation is attributed to the presence of events 
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Figure 3. Final charged particles MD for the 2-jet 
sample of events (y m i n = 0.02) at y/s = 91 GeV 
from the Delphi Collaboration (diamonds) and 
fits with eq. |7 (dotted line, open circles in the 
residuals) and 10| (solid line, filled circles in the 
residuals). 



with quarks of heavy flavour. This suggests that 
longitudinal observables can reveal the flavour 
structure of the original hard event. Second, the 
result by Delphi on the MD in one hemisphere 
(j)) : a sample enriched in bb events was found to be 
essentially identical in shape to the MD of the full 
sample, apart from a shift of one unit. If taken 
literally, and if forward-backward correlations are 
neglected except for the requirement of even total 
charged multiplicity, this result suggests that the 
bb events MD, when the fraction of bb events is 
small, is very similar to the light qq MD, except 
for a translation of two units. 
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Figure 4. H q ratio from the MD's in figure || 



Guided by all the above considerations, we pro- 
pose to use a NBD to describe the MD of a 2-jet 
event of fixed flavour. We identify two types of 
hadronic events: those originating from a bb pair 
and those from a qq pair, where with q we indi- 
cate all flavours lighter than b. We thus use for 
the MD of 2-jet events the weighted sum of two 
NBD's: 

Pfi ^hcavy^n ^(/^heavy? ~t~ 

(l-a hcavy )Pi NBD) (Wight,£0 (10) 

Here «hcavy = 0.22 is the experimentally deter- 
mined fraction of bb decays of the Z° at s/s = 91 
GeV. The subscripts 'heavy' and 'light' refer re- 
spectively to bb and qq events. It should be no- 
ticed that the parameter k is the same in both 
NBD's, while the difference between the average 
values is not fixed. This reflects the above con- 
siderations on single hemisphere multiplicity. 

The description of published data by the above 
parameterization is very good, as shown in table 
0d and in figures |^ and || (solid lines and solid 
dots), from all points of view: chi-square, resid- 
uals and H q analysis. It is further supported by 
the fact that a few simpler alternatives do not 
work as well |l^| . 
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4. SUMMARY AND CONCLUSIONS 

A relatively simple parameterization, based 
on phenomenological considerations and previ- 
ous experimental analyses and results, has been 
shown to describe very well the charged particle 
multiplicity distribution in e + e~ annihilation at 
LEP energy, including the tail, that is the mo- 
ments of high order. This is so even at the very 
specific level of the 2-jet sample of events. 

This parameterization is characterized by the 
presence of two types of events (heavy and light 
quarks decays of the Z°) and by the fact that each 
of these is described by a negative binomial dis- 
tribution. Since for 2-jet events each hemisphere 
correspond to one jet, and the two hemisphere 
are practically uncorrelated, we conclude that the 
MD of a "single quark-jet" is described by the 
NBD. Thus we are lead to suggest that the NBD 
is the fundamental building block of multiparticlc 
production in e + e~~ annihilation. 

Finally, it was shown that the single quark-jet 
NBD has a parameter k which is the same for all 
flavours. Because of the relation between k and 
K 2 , eq. H above, we conclude that true correla- 
tions are expected to be flavour independent. 

It should be stressed that all mentioned con- 
clusions are predictions which can easily be the 
subject of experimental testing; we think indeed 
that the results of these tests will improve our un- 
derstanding of multiparticlc dynamics and of the 
properties of strong interactions. 
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